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9 Hi Abstract 
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^ . We begin with the discussion of the relationship between the trace of the energy 

momentum tensor and the gluon condensate at finite temperatures. Using the recent 
numerical data from the simulations of lattice gauge theory for quantum 
chromodynamics(QCD) we present the computational evaluations for the gluon 
condensate. A short discussion of the properties of deconfinement and the implications 
on the high temperature limit are included. We also mention the case of the massive 

quarks where some of the properties of the condensate appear to change. We put 
together these results with some ideas related to the dilatation current. We draw the 
conclusion that the nature of the strong interactions implies that the thermodynamics of 
quarks and gluons never approach even at very high temperatures that of an ideal 

ultrarelativistic gas. 
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1 Introduction 

Out of the present findings based on previous knowledge new science is revealed. 



This statement particularly pertains to the recent experimental results on the Bose- 
Einstein condensation of atomic alkalai vapors [1]. Although this talk will not be directly 
related to these results, it draws much inspiration from them. 

In this presentation we discuss the consequences for the gluon condensate at finite 
temperature of the recent high precision lattice results [2, 3] for the equation of state in 
lattice gauge theory. We compare these results with various other calculations of the 
expected high temperature behavior. The essential relationship is the trace anomaly that 
is due to the scale variance of quantum chromodynamics (QCD). It relates the trace of 
the energy momentum tensor to the square of the gluon field strenths through the renor- 
malization group beta function. Here we shall expand upon the approach investigated 
in [4], for which the consequences of the new finite temperature lattice data for SU(N C ) 
gauge theory for the gluon condensate [5] have been presented. 



2 The Trace Anomaly at finite Temperature 

The idea of the relationship between the trace of the energy momentum tensor and the 
gluon condensate has been studied for finite temperature by Leutwyler [6] in relation to the 
problems of deconfinement and chiral symmetry. He starts with a detailed discussion of the 
trace anomaly based on the interaction between Goldstone bosons in chiral perturbation 
theory. Central to his discussion is the role of the energy momentum tensor, whose trace 
is directly related to the gluon field strength. It is important to note that the energy 
momentum tensor T^ U (T) can be separated into the zero temperature or confined part, 
Ttf u , and the finite temperature contribution # Miy (T) as follows: 



The zero temperature part, Tq V , has the standard problems with infinities of any ground 
state. It has been discussed by Shifman, Vainshtein and Zakharov [7] in relation to the 
nonperturbative effects in QCD and the operator product expansion. In what follows we 
shall just use a bag type of model [8] as a means of stepping around these difficulties 
since we are only interested here in the thermal properties. The finite temperature part, 
which is zero at T = 0, is free of such problems. We shall see in the next section how the 
diagonal elements of 6^ V (T) are calculated in a straightforward way on the lattice. The 
trace 9^(T) is connected to the thermodynamical contribution to the energy density e(T) 
and pressure p(T) for relativistic fields [9] and relativistic hydrodynamics [10] 



(2.1) 



0£(r) = c(r)-3p(r). 



(2.2) 
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The gluon field strength tensor is denoted by G% v , where a is the color index for SU(N). 
The basic equation for the relationship between the gluon condensate and the trace of the 
energy momentum tensor at finite temperature was written down by Leutwyler [6] using 
the trace anomaly in the form 

(G 2 ) T = (G 2 ) - <0>, (2.3) 
where the gluon field strength squared summed over the colors is 

G 2 = ^i^G^G* (2.4) 



for which the brackets with the subscript T mean thermal average. The renormalization 
group beta function (3(g) in terms of the coupling may be written as 

m = = -7^ {UNc - 2Nf)g3 + 0{g5) - (2 - 5) 

He has calculated [6] for two massless quarks using the low temperature chiral expansion 
the trace of the energy momentum tensor at finite temperature in the following form: 

«5K = -mii{ ln Y ~ i} + 0(rIO) - < 2 ' 6) 

where the logarithmic scale factor A p is about 0.275GeV and the pion decay constant 
F n has the value of 0.093GeV . The value of the gluon condensate for the vacuum (G 2 )o 
was taken to be about 2GeV/ fm 3 , which is consistent with the previously calculated 
values [7]. The results sketched by Leutwyler at QUARK MATTER'96 in Heidelberg [6] 
show a long flat region for (G 2 )t as a function of the temperature until it arrives at values 
of at least 0.1 GeV where it begins to show a falloff from the vacuum value proportional 
to the power T 8 . 



3 Lattice Data for the Equation of State 

The lattice calculation at finite temperature proceeds (talk by Krzysztof Redlich [11]) in 
the following way. From the action expectation value at zero temperature, Po, as well 
as the spatial and temporal action expectation values at finite temperature, P a and P T 
respectively and N T the number of temporal steps, the dimensionless interaction measure 
A(T) [12-14] is given by 

dcr 2 

A(T) = -QNN^a^L- [2P - (P a + P T )\ . (3.1) 

The crucial part of these recent calculations is the use of the full lattice beta function, 
/3f n = adg~ 2 /da in obtaining the lattice spacing a, or scale of the simulation, from the 
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coupling g 2 . Without this accurate information on the temperature scale in lattice units 
it would not be possible to make any claims about the behavior of the gluon condensate. 
The dimensionless interaction measure is equal to the thermal ensemble expectation value 
of (e — 3p)/T 4 . Thus by the equation (2.2) above is equal to the expectation value of the 
trace of the temperature dependent part of the energy momentum tensor, [4] which may 
be written (after suppressing the brakets)as follows: 

0£(T) = A(T) x T 4 . (3.2) 

There are no other contributions to the trace for QCD on the lattice. The heat conduc- 
tivity is zero. Since there are no non-zero conserved quantum numbers and, as well, no 
velocity gradient in the lattice computations, hence no contributions from the viscosity 
terms appear. For a scale invariant system, such as a gas of free massless particles, the 
trace of the energy momentum tensor, equation (3.2), is zero. A system that is scale vari- 
ant, perhaps from a particle mass, has a finite trace, with the value of the trace measuring 
the magnitude of scale breaking. At zero temperature it has been well understood from 
Shifman et al. [7] how in the QCD vacuum the trace of the energy momentum tensor 
relates to the gluon field strength squared, G$. Since the scale breaking in QCD occurs 
explicitly at all orders in a loop expansion, the thermal average of the trace of the energy 
momentum tensor should not go to zero above the deconfinement transition. So a finite 
temperature gluon condensate G 2 (T) related to the degree of scale breaking at all temper- 
atures, can be defined to be equal to the trace. We have used [5] the lattice simulations 
[2,3] in order to get the temperature dependent part of the trace and, thereby, the value 
of the condensate at finite temperature. In what follows we will use the assumptions of 
a bag type of model [8] which includes a bag pressure —B and an energy density of the 
confined state of the same magnitude B. These assumptions mean that this contribution 
to the trace of the energy momentum tensor becomes simply = 4B. Using this form 
we may write the equation of state for the total energy density e(T) and pressure p(T) as 
follows: 

e(T) - 3p(T) = T;(T) (3.3) 

The trace of the energy momentum tensor as a function of the temperature is shown in 
Figure 1. We notice that for T < T c it remains constant at AB. However, above T c in 
both cases there is a rapid rise in Tjf(T). Accordingly, the vacuum gluon condensate Gg 
becomes just 4B. Here we have assumed the value 4B = 0.012GeV A for both cases [7]. 
By taking the published data [2,3] for A(T), and using equations (2.3) and (3.2) we 
obtained the gluon condensate G 2 {T) as shown in Figure 1. 



4 High Temperature Behavior of the Equation of State 

In this section we present a discussion of the properties of the equation of state at tem- 
peratures above 2T C up to around 5T C . If one were to consider a simple bag type of 
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model with the bag constant B independent of the temperature, one finds immediately 
the simple equality 

e - 3p = AB. (4.1) 

This form is in clear conflict with what we found above as shown in Figure 1(a), for which 
there was a steady rise in in its value. Some of the early work in this temperature range was 
done using perturbative estimates by Kallman [15] and Gorenstein and Mogilevsky [16]. 
These authors found essentially a linear rise in the equation of state a function of the 
temperature. Also Montvay and Pietarinen [17] looked at the asymptotic properties of 
the gluon gas. The asymptotic behavior of the form of A(T) assumed by Kallman [15] 
and Gorenstein and Mogilevsky [16] are qualitatively similar to our plots in as far as the 
data extends. However, the latter [16] represent their data in terms of energy density 
and the pressure separately showing how the two curves converge to eachother at small 
values of 1/T 3 . All these results are quite similar to the ideal relativistic massive gas. 
The equation of state is then 

e - 3p = l^TK^). (4.2) 

We can readily see form the asymptotic properties of the modified Bessel function K v (x) 
in this equation that the equation of state is always positive and growing as a function of 
the temperature. In the high temperature limit the equation of state relates to (mT) 2 . 




Figure 1: The left plots show T£(T) for the lattice gauge theories SU(2) and SU(3) as 
indicated. The right plots show the corresponding gluon condensates. The values for the 
critical temperature T c = 0.290,0.264 GeV for SU(2) and 577(3), respectively. 
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We plot here for the sake of comparison with these authors the high temperature 
behavior of A(T) at temperatures in the above range. The following two plots show the 
high temperature behavior of SU(2) and 577(3) explicitly as a function of 1/T 3 . We see 
from these plots that the general appearance is much the same in all the different cases. 
Roughly speaking we could interpret A(T) as a linearly increasing function in the variable 
1/T 3 . In Figure 2(b) we have provided three lattice sizes for the sake of comparison. We 
notice that there exists a considerable gap with no indication of A(T) approaching zero 
at small 1/T 3 . Although we are not able to exactly determine the asymptotic form of the 
function at high temperature, we are able to conclude about certain behavior from given 
models in the range of temperature between 2T C and 5T C . 

Finally we briefly mention the results of finite temperature perturbation theory [18]. 
From the expression for the free energy one can derive the interaction measure in the form 

A = ^C\d*c?(± + i(^) V2 + 0(a)), (4.3) 

where a is related to the coupling by g 2 /47r. The constants Ca and cIa relate to the 
number of colors as N c and — 1 respectively. The significance of the factor y was 
pointed out to us by Bo Andersson [19] as the value in the rapidity of the gluon splitting. 
Thus we see from (4.3) that the temperature dependence of the equation of state is very 
strong in comparison with our plots. 




Figure 2: Figure (a) shows A(T) as a function of 1/T 3 for SU{2) lattice gauge theory 
[2], where the scale of the abscissa is T c ~ 3 x 10~ 3 ; Figure (b) shows A(T) for S77(3) with 
the three lattice sizes [3] 16 3 x 4 (top), 32 3 x 6 (middle), 32 3 x 8 (bottom). The scale of 
the abscissa is the same as for SU(2). 
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5 Gluon Condensate with Quarks 



In this section we would like to discuss the changes due to the presence of dynamical 
quarks with a finite mass. There have been recent computations of the thermodynamical 
quantities in full QCD with two flavors of staggered quarks [20,21], and with four fla- 
vors [22] . These calculations are not yet as accurate as those in pure gauge theory for two 
reasons. The first is the prohibitive cost of obtaining statistics similar to those obtained 
for pure QCD. So the error on the interaction measure is considerably larger. The second 
reason, perhaps more serious, lies in the effect of the quark masses currently simulated. 
They are relatively heavy, which increases the contribution of the quark condensate term 
to the interaction measure, A m (T) in the equation (5.3) below. 

In the presence of massive quarks the trace of the energy-momentum tensor takes the 
form from the trace anomaly [24] as follows: 

(eg) = mSM + (g 2 ), (5.i) 

where m q is the light (renormalized) quark mass and ip q , ip q represent the quark and 
antiquark fields respectively. For the sake of simplicity we choose two light quarks of the 
same mass m q = 6MeV and respecting isospin symmetry. We are now able to write down 
an equation for the temperature dependence of the gluon condensate including the effects 
of the light quarks in the trace anomaly in the following form: 

(G 2 ) T = (G 2 ) + m q $Mo - rnSMr ~ {%)t- (5.2) 

It is possible to see that at very low temperatures the additional contribution to the 
temperature dependence from the quarks is rather insignificant. However, in the range 
where the chiral symmetry is being restored there is an additional effect from the term 
{ip q ip q )T- Well above T c after the chiral symmetry has been completely restored the 
only remaining effect of the quark condensate is that of the vacuum. This term would 
contribute negatively to the gluon condensate. Thus we expect [5] that for the light 
quarks the temperature dependence can be quite important around T c . In fact we may 
rewrite the above equation using the known vacuum estimates for (G 2 )o and m q (ijjip) in 
the simpler form 

(G 2 ) T = B A - m q (^ q ) T - A m (T)T 4 , (5.3) 

where the constant B 4 has the value 0.01183GeV 4 for the u and d quarks alone [7]. The 
use of this value should not imply the given accuracy- in reality the accuracy is still the 
known value of 0.012GeV 4 . However, it should imply that it could be possible directly in 
the neighborhood of T c that the temperature dependence from the quarks is sufficiently 
strong so that the rise in this contribution from temperature compensates to a certain 
extent for the strongly negative tendency of the pure gluon condensate [5]. For the case 
of pure SU(N) we know that the values just below T c of A(T) x T 4 are very small-that 
is, about the same size as m q {ip q ip q )o. Thus, to the extent that the chiral symmetry has 
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not been completely restored, its effect on (G 2 )t will be present below T c . 

As an end to this sketch of the gluon condensate in QCD we will mention a few 
obvious points. Where in simulations on pure SU(N C ) gauge theory we could depend on 
considerable precision in the determination of T c and A(T) as well as numerous other 
thermodynamical functions, this is not the case for the theory with dynamical quarks. 
The statistics for the measurements are generally smaller. The determination of the 
temperature scale is thereby hindered so that it is harder to clearly specify a given quantity 
in terms of T. Thus, in general, we may state that the accuracy for the full QCD is way 
down when compared to the computations of the pure lattice gauge theories. However, 
there is a point that arises from the effect that the temperatures in full QCD are generally 
lower, so that A(T) x T 4 is much smaller [5]. Here we can only speculate with the present 
computations [20-22]. Nevertheless, there could be an indication of how the stability of 
the full QCD keeps (G 2 ) T positive for T > T c . The condensates in full QCD have 
also been considered by Koch and Brown [25]. However, the lattice measurements they 
used were not obtained using a non-perturbative method, nor was the temperature scale 
obtained from the full non-perturbative beta-function. 



6 Dilatation Current 

We have noticed in the previous sections that the fact that the trace of the energy mo- 
mentum tensor does not vanish for the strong interactions has important implications for 
the equation of state. Here we shall discuss briefly some more theoretical results relating 
to T£. This arises with respect to the dilatation current as well as the special conformal 
currents, which are not conserved. The dilatation current may be defined in terms of 
the position four- vector x M and the energy momentum tensor T^ v as simply rr a T MQ . In 
the case of general energy momentum conservation one can find [26] 



= (6.1) 

We now look into a volume in four dimensional space-time V4 containing all the quarks 
and gluons at a fixed temperature T in equilibrium [27]. The above equation (6.1) holds 
when the energy momentum and all the (color) currents are conserved over the surface 
8Va of the properly oriented four- volume V4. Then from (6.1) 

/ V^dS" = [ T»dV 4 , (6.2) 

where the dyxle three-form is T>^dS^ on the three dimensional surface dV^ dual to D^dx^ 
the dilatation current one-form [28]. On the right hand side of (6.2) the integrated form 
J Vi T£dV± is an action integral involving the equation of state. Since T£ > 0, the action 
integral is not zero. This action integral gets quantized with the fields. 
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7 Summary and Conclusions 



Absolutely central to our development of this discussion is the nontrivial equation of 
state (3.3) with Tg(T) > for all reachable temperatures. Thus a strongly interacting 
gas of quarks and gluons in equilibrium is never an ideal gas. 

The lattice simulations for pure SU(N C ) provide not only an accurate computation of 
the interaction measure A(T), but also of the temperature scale due to the calculation 
of the beta function at finite temperature. These properties are all needed in the com- 
putation of 9fc(T) in the previous sections. It is clear that the gluon condensate becomes 
negative in pure gauge theory, and keeps dropping with increasing temperature. We have 
also given a discussion of these quantities in relation to full QCD. In all cases considered 
above we may conclude T^(T) always remains positive so that the divergence of both 
the dilitation and the conformal currents remain positive, and are thereby not conserved. 
Therefore, the scale and conformal symmetries remain broken at all temperatures for 
nonabelian lattice gauge theories after quantization. 

An important point for the physical meaning of these results is the comparison of the 
form A(T) x T 4 with the corresponding term in Bjorken's equation of state pAi [29]. The 
definition of Ai relates to the change of the effective number of degrees of freedom for 
which a first order phase transition may appear. We also notice the general similarity 
in shape of A 1 to the lattice interaction measure with the main difference being that 
the peak in the lattice simulations [2, 3] is slightly above T c . Further comparison with 
Bjorken's model could lead to a better physical understanding of the equation of state 
gotten from lattice simulations. 
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